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Abstract 

Full shell-model diagonalization has been performed to study the structure of neutron-rich nuclei 
around 20 C. We investigate in detail the roles played by the different monopole components of the 
effective interaction in the evolution of the N - 14 shell in C, N and O isotopes. It is found that 
the relevant neutron-neutron monopole terms, V^" jds 2 and V"" 2Sip , contribute significantly to the 
reduction of the N - 14 shell gap in C and N isotopes in comparison with that in O isotopes. The 
origin of this unexpectedly large effect, which is comparable with (sometimes even larger than) 
that caused by the proton-neutron interaction, is related to the enhanced configuration mixing in 
those nuclei due to many-body correlations. Such a scheme is also supported by the large B(E2) 
value in the nucleus 20 C which has been measured recently. 
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1. Introduction 

The study of neutron-rich nuclei with unusually large N/Z ratios is challenging the conven- 
tional view of nuclear structure IHI2]. It is established that, when going from the /^-stability 
line to drip line, the shell structure evolves and new magic numbers may emerge due to the dy- 
namic effects of the nucleon-nucleon interaction. The influence of the proton-neutron interaction 
and its higher order term (namely the tensor force) has been extensively analyzed in recent pub- 
lications IH |3r(^|. One may expect that other components of the nucleon-nucleon interaction 
contribute also to the shell evolution. However, this issue has not been much covered, which 
may be related to the fact that the isoscalar channel of the effective interaction is usually much 
stronger than the isovector part. In Ref. H, it was found that a modification for the isovector 
T — 1 channel of the monopole interaction is needed to reproduce the energies and electromag- 
netic transition properties of low-lying states in nuclei around 17 C. 

In this paper we analyze in detail contributions from the different terms of monopole inter- 
action to the evolution of shell structure. In particular, we will show that the neutron-neutron 
monopole interaction can have significant effect on the evolution of N - 14 shell in carbon iso- 
topes. The rapid decrease of the E{2\) in N = 14 isotones indicates that the N - 14 shell gap 
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existing in oxygen isotopes erodes in nitrogen and carbon isotopes 

JllEMl- At AT = 14, the 

B(E2) value measured through the lifetime of the 2\ state in 20 C B15I1 is much larger than the 
values in 1618 C |T^, 17] and 22 |1_8] while an inelastic scattering measurement presents a much 
smaller B(E2) value in 20 C lfl9ll . In this paper, we will also analyze E2 transition properties in 
neutron-rich carbon isotopes from the viewpoint of collectivity. 

In Sec.|2]we briefly introduce the theoretical framework. The influence of different monopole 
terms of the effective interaction on the evolution of the N - 14 shell is discussed in Sec. 13. II In 
Sec. l3.2l we analyze E2 decay properties of the nuclei I6 . 18 ' 20 Q The calculations are summarized 
in Sec. |4] 



2. Theoretical framework 

The model space we choose contains p and sd shells lEoll . To facilitate the calculation, we 
restrict the maximum number of two for nucleons that can be excited from p to sd shell. This 
truncation is denoted as 2hu> space as usual. It should be mentioned that the full 2hco space 
should also include contributions such as the excitation of one nucleon from the p shell to the fp 
shell or from the Qs shell to the sd shell. Such configurations are neglected in the present study, 
since they do not show any significant influence on the structure of nuclei of concern. Similarly, 
in the present work, the Oftu> truncation implies that no nucleon is excited from the p shell to 
the sd shell. The well-established WBT [21] and MX H22I1 effective Hamiltonians are used. The 
WBT interaction is constructed based on the USD [23] interaction. The p-shell and psd-cross- 
shell matrices in WBT are determined by fitting experimental data. The sd, p, psd and ppsdsd 



parts in the MK interaction are taken from the PW interaction [24], CK interaction |25|], MK 



potential 12211 and Kuo-Brown G matrix 1 26], respectively. Calculations are carried out with a 
newly-established parallel shell-model code described in Ref. 027I1 . 

In 2tioj calculations, one has to remove the spurious states which are caused by the center- 
of-mass motion. As suggested in Ref. B28I1 . the shell model Hamiltonian H is modified to be 
H' = H + BH cm where H cm is the center-of-mass Hamiltonian. One can separate the low-lying 
states from the spurious states by taking a large positive B value, which can move the center-of- 
mass excitations to high excitation energies. We take B = 100 MeV in the present work, which 
is used in some previous works in this region II29I-I3 ill . 

The effective single-particle energy (ESPE) 13211 . describes in a simple manner the evolution 
of shell structure within the shell model framework The ESPE of a given / orbit is 

written as HQ 

f 

where s c ° le is the single-particle energy with respect to the core, {iff\Nj> is the occupancy of 
nucleons in the / orbit and V)/ is the monopole interaction [33]. As in Ref. Jit], one has 

V jf , T = J [1 - (-1)^™^/)] {2j f^y + l) mV\jfhr, (2) 

where (jj'\V\jj')jT is the two-body matrix element and J and T are the angular momentum 
and the isospin of the corresponding two-particle state, respectively. To explore the influence of 
configuration mixing on the ESPE, we will adopt the occupancy in Eq. ([]]) obtained from the full 
shell-model diagonalization. 




Neutron Number 



Figure 1: Calculated £(2|) and ESPEs in even-neutron C, N and O isotopes. In N isotopes, the effective E(2\) are 
deduced from the 27+1 weighted average of E(3/2j) and E(5/2j) 11211 . Experimental data are fromRefs. H0Ul2Ul8t 

m. 



3. Calculations and discussions 



3.1. N=14 shell closure 

In the independent particle model, the last neutrons in 16 C, 18 C and 20 C should occupy the 
y0a?5/2 orbit as in O isotopes. However, the ground state and the first excited state of 15 C are l/2 + 
and 5/2 + states, respectively, which indicates that the v\s\/2 and orbits are inverse in 15 C 

with respect to 17 0. As noted in Ref. lf34ll . 16 C shows a strong (vLsi/2) 2 configuration. The rapid 
decrease of E(2\) in /V = 14 isotones from O to C |I3-12] also indicates that shell structures in 
O and C isotopes are different. 

In order to study the structure difference between O and C isotopes, firstly we present £(2|) 
and ESPEs in the even-neutron C, N and O isotopes, as shown in Fig. Q] For N isotopes, the 
energies are effective £(2|) which are deduced from the 27+1 weighted average of E(3/2*) 
and £(5/2|) 11211 . It is seen that the calculated E{2\) values with WBT and MK interactions 
are systematically larger than experimental values in C and N isotopes. It has been pointed out 
that a reduction of V" n in the sd part of WBT can improve the results ifTH fl2Tl . For such an 
improvement, the modified versions of WBT have been suggested by multiplying the sd part of 
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V m in WBT by a factor of 0.75 for C isotopes (named WBT* [11]) and 0.875 for N isotopes 
(named WBTM [12]). We can see that the E{2\) calculations are improved with WBT* for C 
isotopes and WBTM for N isotopes, shown in Fig.Q] From ls O to 16 C, the quenching of V" is 
explained via the core polarization effect S. In Fig. cn we do not show the 2Hoj results, since 
they are very similar to the 0hco calculations. Both the energy and the vQds/2 - vlsi/2 gap in N 
isotopes are almost in the middle of those in C and O isotopes. 

In C isotopes, the calculated energy gap between vis 1/2 and vOd^ ji orbits is much smaller than 
what is shown in Ref. Ill ill . By calculating ESPEs with shell-model occupancy, the correlation 
effect is approximately included. So the ESPEs are different from those calculated with simple 
occupancy like assuming six neutrons all occupying v0o?5/2 orbits in 20 C and 22 0. The WBT* 
and WBTM give smaller vOc/g/2 - vlsi/2 gaps for C and N isotopes than those by WBT. 

The role that each interaction plays in evolution of N - 14 shell can be investigated based on 
present ESPEs. The energy gap between vlsi/2 and vOds/2 is written as 

Aejv=i4 = s v ij 1/2 - £vOd 5 /2- (3) 
The difference between 2l) C and 22 is written as 

A ejv=14 ( 22 O- 20 C) = A^ =14 ( 22 O)-A e)v=1 4( 20 C) 

= [£vl. Sl/2 ( 22 0) - £vl, I/2 ( 2 °C)] - [£vfl rf5/2 ( 22 0) - £v0d 5/2 ( 2 °C)] 

= 2^[<vlii/2j"|V|vlii/2/)/r - {v0d 5/2 j'\V\v0d 5 /2j')jT] X 

f 

[N r e 2 o)-N r e°c)]. (4) 

The first part is the monopole interaction between the j'-th orbit and vL?i/2 (vOo^) orbit, 
while the second part is the difference of occupancy on j'-th orbit between 22 and 2() C. Besides 
two Op 1/2 protons, the numbers of protons and neutrons on other orbits are also different between 
22 and 20 C. The contribution from each term of the monopole interaction can be analyzed 
by Eq. (0). There are totally 20 terms of the monopole interaction which contribute in 2ha> 
calculation. In Qtiu calculation, the number of the terms reduces to 10 because no protons are in 
sd shell and neutrons in p shell are fully occupied in 22 and 20 C. 

Fig. [2] presents AeAi = i 4 ( 22 0- 21 N) and Ae/ V= i 4 ( 22 O- 20 C), respectively. It is seen that the 
proton-neutron interaction, V^'" (/ - Vp" 2Sl/2 , has large contribution to Ae,v = i4 as discussed in 
Ref. [11]. The two terms, V™ d ^ and V^" , have also significant contributions to Asn=u- In 
2Hu> calculation, the influence of other 16 terms is small. When removing two protons from 22 0, 
the proton-neutron interaction reduces the N = 14- gap and enlarges the mixing between vlsi/2 
and vOt/5/2 orbits. At the same time, the neutron-neutron interaction contributes to the shell 
evolution because the neutrons in vlsi/2 and vQds/2 orbits rearrange as a result of correlation 
effect from shell-model diagonalization. Although the ESPEs do not entirely include the corre- 
lation effect because the multipole interaction is not taken into account, we can approximately 
study the role of neutron-neutron correlation in shell evolution by applying occupancy obtained 
from full shell-model diagonalization to Eq. (Q]l. The attractive neutron-neutron interactions are 
mostly due to the attractive singlet-even (SE) channel of the central force fl. Specially, V"" iSi 2 
is completely from this channel. 

The study of shell evolution extends our understanding of nucleon-nucleon interaction, espe- 
cially the tensor force [Ql . From a general view of shell evolution, when proton (neutron) 
number changes in isotones (isotopes), the neutron (proton) shell structure evolves due to the 

4 




Figure 2: Contribution of each interaction in N = 14 gap difference between 22 and 21 N (left panel), 22 and 20 C (right 
panel). The dashed line is for the value obtained from observed levels in 22 and 21 N 1121 . The 2Tm> calculations are 
indicated, otherwise QHai calculations. 



proton-neutron interaction. Here we present a detailed investigation on N - 14 shell evolu- 
tion to point out that the role played by the proton-neutron interaction in shell evolution can 
be "enhanced" by the attractive neutron-neutron interaction due to the many-body correlations. 
It is important to notice this "enhancement" if one wants to study the nuclear force from shell 
evolution. It is interesting to investigate the contribution of V p ? and V"" in other cases of shell 
evolution in order to gain a complete understanding of how nuclear force drives the shell evolu- 
tion. The N - 16 shell existed in O isotopes will dramatically shrink in Si isotopes [3]. We adopt 
the similar method which we use in the study of N - 14 shell for a detailed investigation of the 
N - 16 shell evolution from 24 to 30 Si. The result shows that, besides the proton-neutron in- 
teraction which was discussed in Ref. [3D, the neutron-neutron interactions, VT j and V"" „ , 

"3/2"3/2 $1/2*1/2' 

contribute to the N - 16 shell evolution from 24 to 3() Si. Because of the limitation of calcu- 
lation, it is hard to perform full shell-model diagonalization to study shell evolution in heavier 
regions, such as that in N = 5 1 isotones which was discussed in Refs. D4j,|5[] . One may expect that 
more information from shell evolution can be gained through full shell-model diagonalization in 
heavier regions in the future. 

3.2. E2 decay properties 

Existing shell-model interactions reproduce well the observed B(E2) values for carbon iso- 
topes II15[ 13711 including the large B(E2) value in 20 C. In shell model, the B(E2) is calculated 
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Table 1: Calculated E2 transition matrix elements and B(E2) values for even-even C. The SFO and SFO-tls calcu- 
lations 1 37] are given for comparison. 
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by 



B(E2, ji -» j f ) = — !— [efA p + e e n ff A n f, (5) 

where /; and are the angular moments of the initial and final states of the E2 transition. A p and 
A n are proton and neutron E2 transition matrix elements, respectively. In TableQ] we present A p , 
A n and B(E2) values for even-even 16 ~ 20 C. In the calculations, the oscillator parameter of hco — 
45A _1 / 3 -25A~ 2 / 3 is employed. An approximate 1/A-dependent effective charge is used l39ll . For 
16 " 2() C,the effective charges are e e / f (e e / f ) = 1.16(0.33), 1.11(0.27), 1.07(0.22), respectively Q. 
It is seen that the A p increases significantly from 16 C to 20 C. Its contribution in B(E2) values is 
multiplied by e p ^ . Thus the B(E2) value in 20 C is much larger than those in 16 C and 18 C. As 
discussed in Refs. |l5, 37], the reduction in the nOpi/2 - tt0/?3/2 gap from 16 C to 20 C would be 
one reason why the proton excitation probability and the B(E2) value increase from 16 C to 20 C. 

The large B(E2) value in 2l) C can also be explained from the view of collectivity. In our 
previous work, we claimed that 2l) C may have large quadrupole deformation 13811 . The intrinsic 
quadrupole moment Qo(s) of the 2| state can be deduced from the spectroscopic quadrupole 
moment Q spec - If the 2| state is a rotational deformed state, the intrinsic quadrupole moment 
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Table 2: Calculated intrinsic quadrupole moments Qq(s) (from spectroscopic property) and |Qo(OI (from E2 transition), 

and the ratio Qo(s)/Qo(t)\. The experimental |Qo(OI are deduced from the observed B(E2) values fl5Hl7ll . 

g^MC/gj \Qo(t)\(fm z ) \Qo(s)/Q (t)\ 

MK-Oftw 10.07 16.63 0.61 

WBT-Oftw 13.25 14.48 0.91 

WBT-2ftw 13.93 13.96 1.00 

WBT*-0/zw 13.53 14.02 0.97 
Expt. 11.43 + 1.98 [17] 

14.44+ 1.27 [16] 

18 C 

MK-OSw 11.02 17.52 0.63 

WBT-OSw 11.19 15.77 0.71 

WBT-2fcw 8.72 15.81 0.55 

WBT*-0fcw 9.21 15.75 0.58 
Expt. 14.70 + 2.05 [17] 

20 C 

MK-OSw -20.14 20.59 0.98 

WBT-Oftw -18.36 18.44 1.00 

WBT-2fcw -17.55 18.05 0.97 

WBT*-0/zw -17.44 17.92 0.97 

Expt. 19.41+^jj [15] 



(written as Qo(t)) can also be deduced from the B(E2; 2| — > 0|) value of the E2 transition. In the 
case of a good deformed rotor, these two definitions of the intrinsic quadrupole moments should 
give the similar values. One has 



(J +1)(2/ + 3) 
3K 2 -J(J+ 1)- 



GoO) = \^ 7/7 , ,[ Qspec{J), (K * 1), (6) 



and 



B(E2;J -» J -2) = — e 2 \ < JK20\(J - 2)K > \ 2 Q (t) 2 , 
Ion 

(K + 1/2, 1). (7) 

In our case, we have K — and 7 = 2. The values of B(E2) and Q spec can be calculated in 
shell model. Table|2]lists the intrinsic quadrupole moments for 16 C, 18 C and 20 C, and the ratios 
between Qo(s) and Qo(t). 

It is seen that Qo{s) and Qo{t) are almost the same in 20 C, which, associated with a negative 
quadrupole moment, would indicate that 20 C may have an oblate deformed structure, and thus the 
large B(E2) value in 20 C might be understood with a consideration of the deformed collectivity. 
It may be of interest to understand why such collectivity is found in a nucleus as light as 20 C. The 
six valence neutrons have strong configuration mixing due to the disappearance of the N - \A 
shell, and, due to the reduction of the proton 0p\/2-0py2 gap in 20 C, the probability of the proton 
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excitation to Op 1/2 orbit increases, as discussed above. The strong mixing of configurations in 
both protons and neutrons would induce the enhanced collective property in 20 C. The results from 
AMD lf40h and deformed Skyrme-Hartree-Fock l39tl have also shown an oblate deformation for 
20 C. Compared with 20 C, the ratio between Qq(s) and Qo(f) in 18 C is significantly less than the 
number of 1 in both WBT and MK calculations. Although we can not exclude other collective 
freedoms in 18 C in the present study, the smaller Qu(s) in 18 C compared with 20 C would indicate 
that the quadrupole deformation in 18 C should be smaller than that in 20 C. In case of 16 C, the 
WBT and MK give rather different results. Calculations with the WBT interaction show that the 
nucleus 16 C is well prolate deformed. 

In the shell model where a spherical basis is used usually, the collectivity (or deformation) 
may be discussed in the language of configuration admixtures. In mean-field models, the col- 
lectivity (or deformation) can be discussed qualitatively by the filling of deformed orbits. The 
deformed Skyrme-Hartree-Fock calculation has been performed for 17 C l4lll . giving coexist- 
ing prolate and oblate minima in each potential energy surface of the ground and lowly-excited 
states. It was predicted that the 3/2 + g.s. has the lowest energy minimum at a prolate deformation 
while the first excited 1 /2 + state has the lowest minimum at an oblate shape [41]. The authors 
explained also that the observed hindered Ml transition between the 1 /2 + and 3/2 + states would 
be attributed to the shape difference of the two states 14111 . If we look at the Nilsson diagram 
of single-particle orbits, the l/2[220] orbit of the Oc/5/2 shell has a prolate-driving force and 
5/2[202] has an oblate-driving force, while 3/2[211] has a prolate-oblate competing force. For 
the 3/2 + g.s. in 17 C, the three valence neutrons outside the N = 8 closed shell would have two on 
the l/2[220] orbit and one on 3/2[21 1], which would lead to a prolate deformation for the ground 
state. The first excited l/2 + state at the lowest energy should have two neutrons on 5/2[202] and 
one on l/2[220], which would make an oblate deformation. Back to the even-even C isotopes 
which we are investigating, the similar qualitative analyses of deformations may be made. The 
16,28,20(-; isotopes should have coexisting prolate and oblate shapes, similar to 17 C. The g.s. of 
16 C would favor a prolate shape with two valence neutrons occupying the prolate l/2[220] orbit. 
20 C would favor an oblate deformation with two neutrons on the oblate 5/2[202] orbit, while 18 C 
might be more shape-competitive between prolate and oblate deformations. These qualitative 
analyses are consistent with the shell-model calculations given in Table[2] However, it should be 
noted that the deformations of the carbon isotopes are rather soft, compared to heavy nuclei. 



4. Summary 

In summary, the nuclei around 20 C have been studied with shell model in Qfico and 2hu> 
model spaces. The N - 14 shell closure appears in O isotopes clearly, but vanishes in C iso- 
topes II 1 lL Il2ll . Besides the proton-neutron interaction which was discussed in Ref. Ill ill , the 
neutron-neutron interactions, V™ and V"" 2Sip , play also important role. When the N - 14 
shell gap is reduced, neutrons in the lsi/2 and (W5/2 orbits rearrange and the interactions V^ 2 ^ S2 
and V™ . can have significant effects on the shell evolution. In 2Hu> calculation, other 16 terms 

L M/2-' 1/2 ' 

of monopole interaction have very small contributions to the shell evolution. It would be useful 
to do systematic investigations on the contributions of V p ? and V"" in order to understand the 
nuclear force from shell evolution. 

As a special case of nuclei near the neutron drip line, we have analyzed the B(E2) value for 
20 C in which experimental data have been available. In Refs. Il5l 13711 . it was pointed out that 
the large B(E2) in 20 C is caused by the reduction of the proton Opi/2 - O/J3/2 gap, which leads 
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to a large proton transition matrix element. The present shell-model calculation gives the similar 
result. Usually, a large B(E2) value is connected with collectivity (or deformation). Therefore, 
we have made a further investigation of the problem. Our calculation gives a large quadrupole 
moment (or large B(E2)) for 20 C, indicating an increased collectivity (or deformation). In the 
shell-model language, the enhanced collectivity is caused by strong configuration mixings in 
both valence protons and neutrons. 
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